Equation of State for Neutralino Star as a Form of Cold Dark Matter 
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In order to study the structure of neutralino star and dark galaxy, we consider dynamical in- 
teractions due to boson-exchange in the neutralino matter. Taking into account interactions of 
neutralinos with bosons, we derive the equation of state (EOS) of neutralino stars in terms of the 
relativistic mean field approach. Then we apply the resulting EOS to investigate properties of the 
neutralino star such as its density profile and mass limit. For example, if the neutralino mass is 
around 1 TeV, the Oppenheimer mass limit of the neutralino star is obtained as 6.06 x 10~ Mq, 
and the corresponding radius is about 7.8 mm. Actually, due to an increasing annihilation rate as 
indicated by our calculation, this dense state can never be realized in practice. Our results also 
show that the low density neutralino star may be a possible aggregation of the cold dark matter. 

PACS numbers: 26.60.-|-c,11.30.Pb,95.35+d 



I. INTRODUCTION 

The SNe la and WMAP data P, [2| in observational 
cosmology indicate that our Universe is composed of ap- 
proximately 4% ordinary matter, 22% dark matter and 
74% dark energy. To explain the nature of the cold dark 
matter is one of the most challenging tasks in modern 
physics. Meanwhile, theorists take great efforts to find 
new physics beyond the standard model (SM). The su- 
persymmetry (SJ is believed to be the most possible ex- 
tension of the SM. Within the framework of the mini- 
mal supersymmetric standard model (MSSM), the light- 
est neutralino x5 is a. natural and attractive candidate for 
the cold dark matter [J] . Annihilation of the neutralinos 
provides some detectable signals, such as the high en- 
ergy gamma rays [5|, Q. The authors of Refs. [3, [S, Q 
proposed a model of "neutralino star" as a gravitation- 
ally compressed object of non-baryonic cold dark mat- 
ter. They calculated the diffuse galactic and extragalac- 
tic gamma radiation and claimed that the results are in 
reasonable agreement with observations. They used some 
phenomenological formulae to describe the density distri- 
bution in the neutralino star. Following the literature, we 
are going to study the properties of the neutralino star 
by concerning the dynamical interactions among neu- 
tralinos. The concept of "neutralino star" in this paper 
represents a compact star that only consists of neutrali- 
nos, and it is slightly different from the picture of Refs. 
0,&i. 

The hierarchical structure of dark matter should be in- 
vestigated if we want to know the nature of the dark mat- 
ter. Physical mechanisms have been proposed for the for- 
mation of dark matter objects in literature. Gurevich et 
al. studied how far does the hierarchical structure of dark 
matter go for small scales Q . The scale length of the cre- 
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ated object is determined by the moment when the Jeans 
instability reaches the nonlinear stage [3]. They also 
show that the dark matter particles can be self-trapped 
by gravitational forces and form dark matter stars [1Q|. 
Kolb and Tkachev studied the nonlinear dynamics for 
axions to form boson stars ITl|.[l3l . Matos et al. studied 
the scalar dark matter [13l . Il4l | and examined the pos- 
sibility that galactic halos are coUisionless ensembles of 
scalar "massive compact halo objects" whose mass is in 
the range m < IQ-^M© or SOMq < to < IOQMq [i|]. 
Later we will show that the mass of the neutralino star 
also falls in this range. Neutralinos as stable dark mat- 
ter can also aggregate to form structures, thus, the self- 
annihilation and aggregation processes caused by gravity 
may give rise to a dynamical equilibrium, and the resul- 
tant equilibrium would lead to a formation of the dark 
matter star. 

Relativistic mean field (RMF) approach [1^ has been 
widely used in nuclear physics and particle physics with 
dense medium, and it provides a successful description of 
nuclear matter and finite nuclei, including unstable nu- 
clei. It has also been applied to construct the equation of 
state (EOS) for neutron stars and supernovae. Therefore 
one can believe that it is applicable to our case. Starting 
with an effective Lagrangian of nucleon-meson interac- 
tions, one can determine the EOS of the neutron star by 
using the RMF theory, and then apply it to study neu- 
tron star properties [13, [iM [iMl- The situation of the 
neutralino star is similar, as both nucleons and neutrali- 
nos are stable fermions except neutrons do not annihilate 
each other, but neutralonos do. Thus, it is quite proba- 
ble that the neutralino star stays in an equilibrium of the 
forces caused by the gravitation and the degenerate pres- 
sure and mutual interactions. In this scenario neutrali- 
nos can interact with each other by exchanging neutral 
Higgs bosons and Z boson, just as the mutual interaction 
between nucleons is mediated by mesons. Narain et al. 
studied the compact stars made of fermionic dark matter 
with interactions, which contribute an additional term in 
the EOS determined with the ideal Fermi gas model, as 
a first order approximation [20|. 



In this work, with the neutrahno-boson interaction, we 
use the RMF approach to construct the EOS of the neu- 
trahno star, and then investigate its properties. By solv- 
ing the Oppenheimer-VolkofF (OV) equation, the mass 
density distribution is obtained. For simphcity, we as- 
sume that the constituents of neutrahno star are only 
the lightest neutralinos, which are interacting with each 
other by exchanging four MSSM bosons, and adopt a 
parameter set for the numerical computation, which is 
carried out in terms of the program DarkSUSY 4.1 [2]| . 
The mass of the neutralino star as a function of its cen- 
tral mass density is given and this result is consistent 
with the OV limit calculated with the ideal Fermi gas 
model [221 • The result of the dynamical simulations of 
the neutralino star shows that the neutralino star is a 
possible existent form of the cold dark matter, especially 
in the dark galaxy. 

The paper is organized as follows: In the next sec- 
tion, we use the RMF approach to construct the EOS of 
the neutralino star. In Sec. Ill, we apply the EOS to 
study the properties of neutralino star. The last section 
is devoted to our conclusion and the discussion on dark 
galaxy. 



II. RELATIVISTIC MEAN FIELD THEORY 
FOR NEUTRALINO MATTER 

Within the framework of MSSM, the neutralino can 
interact with four bosons, iJ", /i°, A°, and Z. Here H^ 
and /i° are the two CP-even neutral Higgs bosons, AP is 
the CP-odd neutral Higgs boson, and Z is a vector boson 
in the SM. The effective Lagrangian can be written as 
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-^d^Hd^'H ~ ^mjjH'^ + -d^.hd^'h - -mlh^ 



+ \d,Ad^^A - \m\A^ ~ \f,.F^'' + \mlZ^Z\X) 

where x^ ^^ ^^ A, and Z^ denote the fields of Xi, H^ i 
/i°, A", and Z, respectively, and F^i, = df^Zy — duZ^. 

Starting with this Lagrangian, a set of Euler-Lagrange 
equations can be derived. The Dirac-type equation for 
the neutralino field is given as 



{^r^^. 



gzi^z^,)x = 0, 



(2) 



where m* = m^ -j- qhH + gnh + gA'y^A is the effective 
mass of neutralino. The field equations for TJ", /i", A°, 
and Z are given as 



d,,d^H 



m]jH = 



-9hXX, 
-QhXX, 

'9AX15X, 

d^F'^-' + mlZ'^ = gzXl^'X- 



df.d'^h + mlh 
df.d'^A + m\A 



(3) 
(4) 
(5) 
(6) 



These equations are coupled partial differential equa- 
tions, which are difficult to be solved exactly. To get 



an approximate solution which can properly describe the 
physical picture, we adopt the relativistic mean field ap- 
proximation [i^. In this approach the boson fields are 
treated as classical fields, and the field operators H, 
h, A, and Z^ are replaced by their expectation values. 
Here we consider an approximation that the bosons are 
treated as static infinite matter to obtain simplified equa- 
tions, where derivative terms in the Klein-Gordon equa- 
tions vanish automatically due to the uniformity of the 
matter. The spatial components of the vector field Z^^ 
vanish since there is no vector current in average, i.e., 
(X7*x) = 0- Hence, the equations for the boson fields are 
reduced to 
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At zero temperature, the equations of the boson fields 
can be written as 
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where kp is the Fermi momentum of neutralino. Indeed, 
the temperature in neutralino star is much lower than any 
concerned energy scale, thus it is safe to take the zero- 
temperature approximation. The pseudo-scalar Higgs A'^ 
has no contribution in the RMF approach. These equa- 
tions should be solved self-consistently, then the energy 
density and the pressure can be obtained from the energy 
momentum tensor |16| . The energy density of neutralino 
matter is given by 
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and the pressure is given by 



(15) 
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To obtain a relation between p and p from p — pikp) and 
p — p{kp), we numerically vary the value of kp in a wide 
range, then construct the EOS of the neutralino matter. 
There are several available models in the program 
DarkSUSY ^1 with simpfified MSSM parameterization 
for studies on supersymmetric dark matter, and here we 
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FIG. 1: The EOS of the neutrahno matter, p and p denote 
the energy density and the pressure, respectively. 



FIG. 2: The mass-radius relation of the neutralino star. Ms 
denotes the mass of the sun. 



use the model JE56A_003796. Table [J presents our input 
parameters for the MSSM and Table |lT] shows the pa- 
rameters in the Lagrangian. The EOS of the neutralino 
matter in the energy density range 10~^ ~ 10^'' TeV/fm"^ 
is shown in Fig. [T] 



TABLI 


3 I: Input parameters for the MSSM sector 


Parameter 


Value 



M2 

mA 

tan/3 

ma 

At/nis 
Ab/ms 



1020.2370 
2894.9567 
753.24107 
14.672504 
4001.3887 
-0.57518824 
0.62302097 



TABLE 11: 


The 


parameter 


set for the Lagrangian 


Parameter 






Value 



m^ (GeV) 
m„ (GeV) 
mu (GeV) 
mA (GeV 
mz (GeV) 
an 
9h 

9A 

gz 



1016.79763 

753.36438 

118.760345 

753.24107 

91.187 

0.3131 

-0.02705 

-0.02726i 

0.00143 



III. PROPERTIES OF THE NEUTRALINO 
STAR 

Applying the EOS of the neutralino matter, we calcu- 
late the neutralino star density profile by solving the OV 
equation 



dp 
dr 

where 



GMp 



P 
1 + - 



M{r) 



M 



Attt'"^ p{r')dr' . 



2GM\ 



(17) 



(18) 



For a given central mass density, which is related to the 
Fermi momentum kp, we solve Eq. (|17p from inside to 
outside until the pressure vanishes, then the density pro- 
file and the mass of the neutralino star can be obtained. 

We take the neutralino mass around 1 TeV for example 
to perform the calculations. The mass of the neutralino 
star as a function of its radius is illustrated in Fig. [21 
which shows that the mass limit of the neutralino star is 
6.06 X 10~^Mq (It is less than the mass of the earth.). 
Also, the radius of a neutralino star whose mass just 
takes the mass limit is only 7.8 mm. The mass of the 
neutralino star as a function of its central mass density 
is illustrated in Fig. [3l and the density profile is plotted 
in Fig. H 

By assuming that a compact object of constituents is 
described by the ideal Fermi gas model, the OV limit is 
calculated in Ref. 1221, 



Mov = 2.9924 x IO^Mq ( IL^) ^-V^^ 



- 



(19) 



where m is the mass of the possible constituents (e.g., 
neutron or neutralino) which form the compact star, and 
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FIG. 3: The mass of the neutrahno star as a function of its 
central mass density. Ms denotes the mass of the sun. 



FIG. 5: The difference between the EOS including interac- 
tions (p) and the one without interactions (po) for neutralino 
matter. 
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FIG. 4: The density profile of the neutralino star. The three 
lines corresponds to three different central mass densities. 



g is the degeneracy factor, i.e., g — 2 for Majorana par- 
ticles and 5 = 4 for Dirac particles. By this formula, the 
mass limit of the neutralino star without interactions is 
Mo^J — 6.05 X 1O~^M0, and the limit calculated in our 
approach is consistent with this OV limit. If we set the 
couphngs to be zero, then the EOS reduces to that for 
ideal Fermi gas of neutralinos. Fig. [5] plots the difference 
for the EOS with and without the interactions among 
neutralinos, and it shows that the effect of the concerned 
interactions is rather small, namely, the gravity and the 
degenerate pressure are overwhelmingly important to re- 
sponsible for the formation of neutralino stars. From 
Eq. (fTB]) , we can see that the scalar Higgs bosons H° and 
/i° have negative contributions to the pressure and pro- 
vide an attractive potential (Va), while the vector gauge 
boson Z has a positive contribution to the pressure and 
provides a repulsive potential (Vr). Obviously, in the low 



density region, the interactions can be neglected. As the 
density goes larger, it becomes Va > Vr, and when the 
density is large enough, Vr > Va finally. 



IV. CONCLUSION AND DISCUSSION 

In summary, we have obtained the EOS of neutralino 
matter at zero temperature in the density range from 
10"*' to 10^'' TeV/fm^ and employed the present EOS 
to calculate the neutralino star properties for a specific 
choice of MSSM parameters. We assume that the com- 
pact neutralino star only consists of the neutralino mat- 
ter and we use the RMF approach, which is obvious a 
reasonable approximation in our case. For example, if 
the neutralino mass is around 1 TeV, the Oppenheimer 
mass limit of the neutralino star is 6.06 x 10~^Mq and 
the corresponding radius is 7.8 mm. From the other 
side, there exist two opposite processes, first, the grav- 
ity causes the neutralinos at vicinity close to a core to 
collapse into a small region as the more matter aggre- 
gates, the denser the star becomes, but meanwhile as 
the star becomes denser, the annihilation rates of neu- 
tralinos would rapidly increase and neutralinos annihi- 
late into SM particles, which would escape from the star. 
The process prevents more neutralinos to collapse into a 
core. The final state of the neutralino star would be a 
result of a competition of the two processes and in our 
later work we are going to evaluate the resultant den- 
sities which would determine the size and average mass 
of such neutralino stars. The resultant stars might be a 
candidate for main components of a dark galaxy. 

One is tempted to postulate that the neutralino star 
is a possible existent form of the cold dark matter, espe- 
cially the relatively low density neutralino star in dark 
galaxy. VIRGOHI 21, a dark galaxy in the Virgo cluster, 



has been discovered recently [23[. With behef on vahd- 
ity of supersymmetry, it is reasonable to assume that the 
neutralino matter might be the main content of the dark 
galaxy. The high energy gamma ray as detectable signals 
produced by dark matter is intensely studied in the lit- 
erature [2J, l2a| • The central density of a neutralino star 
near its mass limit is around 10^^ GeV/cm"^, 53 orders 
of magnitude higher than the local dark matter density 
of smooth profile. If the compact neutralino star were a 
main component in the dark galaxy, not only its annihi- 
lation would dramatically change the profile of the dark 
galaxy, but also the gamma ray flux would contradict the 
experimental data. However, we should not exclude the 
possibility of the existence of a few neutralino stars espe- 
cially low density ones. For example, the mean density 
of the dark galaxy can be 0.3 GeV/cm'^ by choosing an 
appropriate central density of such stars, which is reason- 
able to explain the observed gamma ray flux. We think 
that the dark matter can accrete on the neutralino star 
due to its strong gravitational force, thus the dynamical 
equilibrium of its self-annihilation and the dark matter 
accretion will cause such stars to live longer. Therefore, 
the compact neutralino star near its mass limit should 
not exist in the dark galaxy, whereas a few small mass 
neutralino stars are a possible aggregation form of dark 
matter. 



As a further study on the interactions among the neu- 
tralino matter more precisely, we should consider not 
only the lightest neutralino x°j but also the heavier ones 
X234, as well as the charginos X12 ^'^'^ other possible 
particles mixed in the neutralino star, even though they 
may decay into lightest neutralino and SM particles. The 
axion is also an interesting candidate to study [ll|, [I4I . 
Generally, for a mixture of neutralinos with bosons or 
heavier fermions the mass limit would be smaller and a 
mixing with lighter fermions will make the mass limit 
larger. In addition, the similar method can also be ap- 
plied to the case of the "neutrino star" [22], thus, the 
properties of the neutron star, neutralino star and neu- 
trino star all can be studied in a unique theoretical frame- 
work. 
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